y ^ ' 7 

NASA TecKriicai Paper 3557 


Improved Method for Prediction of Attainable 
Wing Leading-Edge Thrust 

Harry W. Carlson, Marcus 0. McElroy, V^endy B. Lessard, and L. Arnold McCullers 



""A 


i 


April 1996 


NASA Technical Paper 3557 



Improved Method for Prediction of Attainable 
Wing Leading-Edge Thrust 


Harry W. Carlson 

Lockheed Engineering & Sciences Company • Hampton, Virginia 

Marcus 0. McElroy and Wendy B. Lessard 
Langley Research Center • Hampton, Virginia 

L. Arnold McCullers 

ViGYAN, Inc. • Hampton, Virginia 


National Aeronautics and Space Administration 
Langley Research Center • Hampton, Virginia 23681-0001 


April 1996 



Available electronically at the following URL address: http.7/techreports.larc.nasa.gov/ltrs/ltrs.html 


Printed copies available from the following: 

NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 
(301) 621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 



Contents 

Abstract 1 

Introduction 1 

Symbols 2 

Present Method Development 3 

Normal Airfoil and Flow Parameter Derivation 3 

Theoretical Two-Dimensional Airfoil Analysis (C p j im = C p vac ) 5 

Equivalent Mach Number Concept (C p ]im ^ C p vac ) 8 

Experimental Two-Dimensional Airfoil Analysis (C p ]„„ Calibration) 9 

Guidelines for Present Method Application to Three-Dimensional Wing Computer Codes 11 

Present Method Validation 13 

Future Limiting Pressure Redefinition With Use of AER02S Computer Code 16 

Concluding Remarks 16 

Appendix A — Strategy for Thrust Factor Data Representation 18 

Appendix B — Navier-Stokes Solutions for Two-Dimensional Airfoils 21 

References 24 

Figures 25 


iii 



Abstract 


Prediction of the loss of wing leading-edge thrust and the accompanying increase 
in drag due to lift, when flow is not completely attached, presents a difficult but com- 
monly encountered problem. A method (now called the previous method ) for the 
prediction of attainable leading-edge thrust and the resultant effect on airplane aero- 
dynamic performance has been in use for more than a decade. Recently, the method 
has been revised to enhance its applicability to current airplane design and evalua- 
tion problems. The improved method (called the present method) provides for a 
greater range of airfoil shapes from very sharp to very blunt leading edges. It is also 
based on a wider range of Reynolds numbers than was available for the previous 
method. The present method, when employed in computer codes for aerodynamic 
analysis, generally results in improved correlation with experimental wing-body 
axial-force data and provides reasonable estimates of the measured drag. 


Introduction 

The loss of wing leading-edge thrust and the accom- 
panying increase in drag due to lift, when flow is not 
completely attached, is a commonly encountered prob- 
lem. Such problems are particularly severe for airplanes 
that employ the thin highly swept wings required to pro- 
vide some degree of supersonic flight capability. Devel- 
opment of a generalized system for the prediction of 
airplane lifting efficiency is dependent on a means of 
evaluating the amount of wing leading-edge thrust that 
can actually be attained. 

Leading-edge thrust results from the high velocities 
and the accompanying low pressures that occur as air 
flows from a stagnation point on the under surface of the 
wing around the leading edge to the upper surface. This 
thrust can be developed at supersonic speeds as well as at 
subsonic speeds provided that the component of the ffee- 
stream Mach number normal to the wing leading edge is 
less than 1.0. The loss of leading-edge thrust resulting 
from flow separation is dependent on flight conditions 
(i.e., Mach number and Reynolds number), on wing sec- 
tion geometry (i.e., maximum thickness, location of max- 
imum thickness, and leading-edge radius), and on the 
incidence of the wing section leading-edge mean camber 
surface relative to the local flow. Flow separation charac- 
teristics of highly swept wings can be especially complex 
because the upwash ahead of the leading edge increases 
dramatically from the wing centerline to regions near the 
wingtip. 

A first step in the development of a practical system 
for the analysis of separated flow problems was provided 
by the theoretical-empirical attainable leading-edge 
thrust prediction method advanced in reference 1. This 
method is based on the following: 

1. Use of simple sweep theory to permit a two- 
dimensional analysis 


2. Employment of theoretical airfoil computer codes to 
define thrust dependence on section geometric charac- 
teristics with pressures limited to a vacuum 

3. Generalization of the thrust dependence on limiting 
pressures to include the more severe limitations of 
realistically achievable pressures 

4. Examination of experimental two-dimensional airfoil 
data to define the more realistic limiting pressure 
dependence on local Mach and Reynolds numbers 

The attainable thrust prediction method was employed in 
computer codes for the estimation of drag due to lift 
(AER02S, refs. 2 and 3) and the design of wing surfaces 
for minimization of drag (WINGDES, ref. 4). The aero- 
dynamic analysis computer code (AER02S) provides 
estimates of the aerodynamic characteristics of a wing 
and body in combination with either a canard or a hori- 
zontal tail surface. It includes provisions for convenient 
handling of deflected wing leading- and trailing-edge 
flaps. This computer code is applicable only for subsonic 
speeds. The wing design computer code (WINGDES) 
provides both a design and an analysis capability and 
is applicable to both subsonic and supersonic speeds. 
However, it cannot accommodate a second surface and is 
not arranged for convenient handling of flaps. Refer- 
ences 5-8 provide an updated description of these com- 
puter codes and examples of their application to practical 
problems. 

Although computer codes employing the attainable 
thrust method have been successfully applied to a variety 
of problems, a need for certain improvements to the basic 
attainable thrust prediction procedures has been recog- 
nized. One deficiency of the method described in refer- 
ence 1 is a poor representation of the thrust-producing 
capabilities of airfoil sections with sharp or nearly sharp 
leading edges. The method, as originally formulated, dic- 
tated that airfoils with a leading-edge radius of zero 
could produce no thrust. As shown by theoretical 



two-dimensional airfoil data, which is discussed later, 
and by experimental investigations (e.g., ref. 9J, sharp 
and nearly sharp airfoil sections do in fact provide appre- 
ciable levels of attainable thrust. In the application of the 
method described in reference 1, the problem was over- 
come in part by substitution of a radius other than zero 
for nominally sharp leading edges; however, this process 
depends on the skill and experience of the computer 
code user. The present method described in this paper 
provides a better solution in which the theoretical 
two-dimensional airfoil matrix is expanded to include a 
leading-edge radius of zero. With this change the method 
is applicable to a continuous range of leading-edge radii 
from zero through the standard values to very large val- 
ues approaching half of the wing maximum thickness. 

Expansion of the two-dimensional airfoil matrix to 
include variations in location of maximum thickness was 
accomplished by a revised relationship between stream- 
wise airfoil sections of the wing and the derived two- 
dimensional sections, a relationship that results in much 
closer representation of the real flow over a lifting sur- 
face. Revision of the attainable thrust prediction method 
also provided an opportunity to take advantage of infor- 
mation relating to the effect of Reynolds number on 
attainable thrust that was not available before publication 
of reference 1. In reference 1, the two-dimensional 
experimental data used to define limiting pressures were 
restricted to R < 8 X 10 6 (based on the chord). The 
present method discussed herein makes use of data 
obtained up to R = 30 x 10 6 . 

Because revisions to the previous method are quite 
extensive, the development of the present method is cov- 
ered in detail, even at the expense of some repetition. 
Some examples of the application of the present method 
to data for wings and wing-body configurations are 
given. Correlations are included for data previously used 
in references 6 and 8 and for new data as well. In addi- 
tion, instructions are given for the evolution of the sys- 
tem to accommodate new two-dimensional airfoil data, 
as it becomes available, so as to provide a more exact and 
more complete formulation of attainable thrust depen- 
dence on Mach and Reynolds numbers. 


C m pitching-moment coefficient 

C N normal-force coefficient 

C p pressure coefficient 

C p jjjj, limiting pressure coefficient used in defini- 

tion of attainable thrust 

2 

C p vac vacuum pressure coefficient, 

y AT 

c local wing chord, in. 

c av average wing chord, ^ , in. 

c A section axial- or chord-force coefficient 

A c A change in section axial- or chord-force coef- 

ficient relative to a = 0° 

c t section theoretical thrust coefficient (from 

linearized theory for zero-thickness airfoils) 

c* section attainable thrust coefficient 

c mean aerodynamic chord, in. 

e i> e 2> e 3 exponents used in curve-fit equation for 

attainable thrust factor 

<? 4 exponent used in curve-fit equation for lim- 

iting pressure coefficient 

K parameter used in curve-fit equation for lim- 

iting pressure coefficient 

K t attainable thrust factor, fraction of theoreti- 

cs c* 

cal thrust actually attainable, — = — — 

C r C t,n 

k parameter used in curve-fit equation for 

attainable thrust factor 

k] ,^ 2 ’^ 3,^4 constants used in airfoil section definition 

M free-stream Mach number 

M e equivalent Mach number replacing M n to 

account for C p>lim * C p>vac 

M n normal Mach number (fig. 2) 


Symbols 


P* 

attainable thrust parameter, 

b 

wing span, in. 



c A 

axial- or chord-force coefficient 



C D 

drag coefficient 

Ptt 

theoretical thrust parameter, 

A C D 

drag coefficient due to lift, C D - C D 0 


<N 

— 
^ 

pj 

Pi 

i-- s T~~ 

C D,0 

drag coefficient at a = 0° for configuration 
with no wing camber or twist 


\cj n yo.5j 
C t,nPn\_ 0.09 


C L lift coefficient 


q dynamic pressure 


2 



R Reynolds number based on mean aero- 

dynamic chord 

r leading-edge radius, in. 

r 

- c n 

r ; - leading-edge radius index, — - 

© 

5 wing area, in. 2 

s spanwise distance, in. 

t section theoretical leading-edge thrust 

t* section attainable leading-edge thrust 

X,Y,Z Cartesian coordinates, positive aft, left, and 
up, respectively (fig. 2) 

x' distance behind wing leading edge 

a angle of attack, deg 



Y ratio of specific heats, 1.4 

8 L n leading-edge flap deflection angle measured 

normal to hinge line, positive with leading 
edge down, deg 

8 r „ trailing-edge flap deflection angle measured 

normal to hinge line, positive with trailing 
edge down, deg 

ri location of maximum thickness as fraction 

of chord 

A le wing leading-edge sweep angle, deg 

t section maximum thickness, in. 

Subscript: 

n quantities pertaining to wing section normal 

to leading edge 

Present Method Development 

Development of the present method for the predic- 
tion of attainable leading-edge thrust is based on the fun- 
damental principle that this force results from pressures 
acting on a surface, and thus, the amount of theoretical 
thrust that can actually be realized is dependent on con- 
straints imposed by wing geometry and achievable pres- 
sure levels. The goal of the present method is to provide 
a means of quantifying that portion of this theoretical 
thrust that can be attained through empirically derived 
mathematical relationships involving wing flow condi- 
tions and wing geometry. The present method, when 
incorporated into existing linearized theory lifting sur- 
face computer codes, estimates attainable leading-edge 


thrust by providing multiplication factors that are applied 
to the calculated theoretical thrust. 

The process used in generation of the prediction 
method is rather complex. Figure 1 provides a visual 
guide to aid in an understanding of the general plan of 
development described in the following paragraphs. 

First, the relationships between wing streamwise air- 
foil sections and sections normal to the leading edge are 
established, and the use of a two-dimensional analysis to 
solve a three-dimensional problem is justified. Simple 
sweep theory is applied to streamwise section flow con- 
ditions, geometry characteristics, and thrust coefficients 
to obtain the corresponding two-dimensional values. 

Second, a computer code for subsonic two- 
dimensional airfoils is employed to define limitations on 
the theoretical thrust imposed by airfoil geometry con- 
straints and by a limiting pressure equal to the vacuum 
pressure. A variety of airfoil sections are treated, and the 
computer code results are collected and represented by a 
relationship expressing the dependence of an attainable 
thrust parameter P* on a theoretical thrust parameter P tt . 

Third, theoretical considerations are employed for 
the general case that enables the vacuum limit thrust 
parameters to be applied to pressure limitations other 
than the vacuum limit. As shown later, this is accom- 
plished by replacement of the normal Mach number with 
an appropriately defined equivalent Mach number. 

Fourth, the actual value of the limiting pressure, 
which reflects the more severe limitations of the real 
flow, is obtained through the use of theoretical and 
experimental two-dimensional airfoil data. For a collec- 
tion of test data, the estimated attainable thrust with 
thrust factors evaluated from the foregoing relationships 
and a trial value of the limiting pressure are compared 
with the experimental thrust data. The calculation is 
repeated with successive trial values of the limiting 
pressure until a sufficiently close match of estimated 
and measured thrusts is achieved. Collected data are 
then empirically represented to establish the variation of 
the limiting pressure with Mach number and Reynolds 
number. 

The complete present method, when applied to lift- 
ing surface aerodynamic computer codes, modifies the 
computed results and replaces the theoretical leading- 
edge thrust with an estimated attainable thrust value. The 
procedure calls for the treatment of a number of span sta- 
tions from the wing root chord to the wingtip. For each 
station, the theoretical thrust is multiplied by a factor 
dependent on the local flow conditions and the local air- 
foil geometry. 
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Normal Airfoil and Flow Parameter Derivation 


The derivation of normal airfoil sections and appro- 
priate flow parameters based on simple sweep theory is 
explained with the aid of figure 2 and sketch A. As 
shown in figure 2, a different sweep angle Aj e may be 
required for each span station. Simple sweep theory is 
applicable to a constant-chord wing of infinite aspect 
ratio as shown in sketch A. 



The fundamental concept is based on constant flow 
properties along lines parallel to the leading edge so that 
the addition of a velocity field with vectors parallel to the 
leading edge will not alter the flow. Thus, the free-stream 
flow, in terms of Mach number, can be divided into the 
following two components: a flow perpendicular to the 
leading edge, which determines the flow properties, and 
a flow parallel to the leading edge, which has no influ- 
ence. This line of reasoning permits flow properties for 
the infinitely swept wing to be calculated from two- 
dimensional flow properties for a section perpendicular 
to the leading edge immersed in a uniform flow of veloc- 
ity M n = M cos A le . 

Even for wings that depart dramatically from a 
constant-chord, constant-sweep condition, the concept is 
useful for the analysis of leading-edge thrust because it is 
dependent primarily on flow behavior in the vicinity of 
the leading edge. As shown by the vectors in figure 2 and 
sketch A, a constant velocity M sin A le added to the nor- 
mal velocity McosAj e gives a resultant flow which 
approximates the actual flow over the wing upper surface 
as represented by the arrows. 

The derivation of normal section geometric charac- 
teristics is explained with the aid of sketch B. In accor- 
dance with simple sweep theory principles, the normal 
section chord is defined as 

C n = C cos A le 



and the normal section thickness is the same as that of 
the streamwise section along lines parallel to the wing 
leading edge. Thus, the normal section thickness ratio is 


gff) 1 

c)„ \cj cos A k 


By inspection, the nondimensionalized location of maxi- 
mum thickness for the normal section is the same as that 
for the streamwise section, i.e., r\ n = q. The leading-edge 
radius of the normal section is derived by an equation 
representing a generalized form of an airfoil thickness 
distribution given in reference 10. A modified form of 
that equation, which was adapted to the purposes of this 
study and which has the normal section symbols of this 
paper, is 



+ k 2 ,n X n 


+ k 


3 , n •* 


' 3/2 

n 


+ k A,n x 


'2 

n 


The first term dominates the shape near the leading edge, 
and the constant kj n defines the leading-edge radius. For 
the normal section near the leading edge 


\ tak l,nJ*n= k l,nJ x ' C0S A ] 


le 


For the streamwise section 


T « kyjx 

So that for x n = T, 

’ 7 cos A i« 

The leading-edge radius (ref. 10) for the normal section 
is given by 

= (*!,/ 

2 


and for the streamwise section by 


r = 


( k 0 2 

2 
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Thus, for the normal section 


The thrust vector in the normal direction is 


_ r n _ r r _n C _ r (^1, n) 1 
c)n ~ c n ~ c r C n ~ c (fei) 2 cos A le 

_ r 1 
C C °S 2 A le 


dt 


n 


dt 

cos A le 


the dynamic pressure in the normal direction is 



M J 


= q cos 2 A le 


These transformations define a wing section that can 
be analyzed for the effect of geometry and pressure 
limitations on the attainable leading-edge thrust by use of 
theoretical and experimental two-dimensional wing 
data. 1 


and the incremental distance in the spanwise direction for 
the normal section is 


cos A le 


The basic premise of the attainable thrust prediction 
method is that an attainable thrust factor K p which is 
derived from the two-dimensional analysis, can be 
applied to the calculated theoretical leading-edge thrust 
for the three-dimensional wing at the corresponding span 
station. Thus, an essential part of the present method 
is a connection between two-dimensional and three- 
dimensional thrust coefficients. The section thrust coeffi- 
cient may be defined as the thrust force per unit dynamic 
pressure, per unit chord, and per unit spanwise distance. 
Thus for the streamwise section, 

- ^£_L 

1 dyqc 

Sketch C aids in describing the derivation of the normal 
section thrust coefficient. 



4n the original development of the attainable thrust prediction 
method (ref. 1), normal section dimensions were defined different- 
ly. That transformation, however, was applicable only to a normal 
section with maximum thickness at the 50-percent chord location. 
The normal airfoil section, as generated in this paper, is more con- 
sistent with simple sweep theory concepts and the actual flow over 
the wing surface. 


The normal section thrust coefficient then becomes 


ds <ln c n 

= dt COS A le 1 

COS A le d y qc n COS 2 A le 

dt J_£_ 1 

” dy V~ cc n cos 2 A le 
c 1 

~ Cf c 2 _ r ~ 

n cos A le 

In the WINGDES and AER02S computer codes, the sec- 
tion thrust coefficient is referenced to an average chord 
c av , which then results in 


c 


t, 


n 


c av 1 

C{ C 2 

n COS Aj e 


The normal section Reynolds number differs from 
the streamwise Reynolds number because of changes in 
the velocity (Mach number) and in the local chord. Thus, 


R=R 


c n^n 

c M 


Theoretical Two-Dimensional Airfoil Analysis 

= £p,v ac) 

For the series of symmetrical two-dimensional air- 
foil sections represented in figure 3, the subsonic airfoil 
program of reference 1 1 was employed to define inviscid 
flow pressure distributions and integrated thrust coeffi- 
cients at a Mach number of 0. The airfoils were defined 
by the following equation: 


t = k 


1, n 


+ x„ + ko .,' 3/2 + lc x 2 

L,n n 3,nn 4, n n 


in which the coefficients were selected to produce 
the required section thickness, location of maximum 
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thickness, and leading-edge radius. Maximum airfoil 
thickness ranged from 3 to 15 percent of the chord. The 
location of maximum thickness varied from 10 to 50 per- 
cent of the chord. A leading-edge radius index 



could be defined as an appropriate measure of the rela- 
tive bluntness of airfoil sections. For the majority of 
NACA airfoils, which span a wide range of thicknesses 
and maximum thickness locations (see ref. 10 for a repre- 
sentative collection), this index varies within a narrow 
range of 0.27 to 0.33. An index of 0.3 was chosen to rep- 
resent a standard leading-edge radius. Other leading- 
edge radius indices of interest included values of 0, 0.15, 
0.60, and 1.2. For a given airfoil, pressure distributions 
and thrust coefficients across a range of normal Mach 
numbers were calculated by application of the Prandtl- 
Glauert rule 

= C ,A = Constant 

to pressure distributions obtained at a Mach number of 0. 
This simple means of handling Mach number effects was 
employed for the sake of consistency with linearized the- 
ory methods used for estimating theoretical thrust. A typ- 
ical pressure distribution for a wing section at lifting 
conditions is shown in sketch D. 


Z 



Values of the integrated section thrust coefficient 
without pressure limitation were calculated from 

c t = \§ C P dz 

The section thrust coefficients were relatively indepen- 
dent of the airfoil thickness and the leading-edge radius 
and were in reasonably good agreement with the two- 
dimensional theoretical value for a zero-thickness airfoil, 
(2tt sin 2 a)/P n , except for a small increase that tended 
to be linearly dependent on the airfoil thickness. The suc- 
tion peak of the theoretical pressure distribution can be 
quite large, often exceeding the vacuum pressure limit 
for a given Mach number. Thus, the theoretical section 
thrust coefficient can be unrealistically high. 

To deter mi ne the effect on the thrust coefficient of 
realistically attainable pressure distributions, the com- 
puter code integration was performed with the pressure 
coefficients limited to values greater than or equal to 
specified pressure coefficients C p i j m , which truncates the 
suction peak as shown in sketch E. This pressure limita- 
tion is intended to account, in an approximate way, for 
two of the factors which limit attainable thrust: the real- 
world failure to attain theoretical peak suction pressures, 
and the tendency of experimental suction peaks to occur 
at a more rearward position on the airfoil. The specified 
limiting pressure coefficieqp were set equal to the vac- 
uum pressure limit -2/y M n . As discussed in reference 1 
and as explained later, repetition of the integration with 


Z 
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more severe real-world pressure limitations is not neces- 
sary. The one set of data with C p ]im = C pwac is sufficient 
because results for other limiting pressures can be 
derived from that data with an appropriate substitution of 
an equivalent Mach number M e for M n . The integrated 
value of c t associated with C p lim is designated c* and 
the attainable thrust factor is simply AT, = c* /c f with 

,o „ r, n r, n 

c t n = (2tt sin a)/p n . A normal section thrust coefficient 
c t n defined in this way is appropriate only for a zero- 
thickness two-dimensional airfoil. The zero-thickness 
reference is chosen because the thrust factor is eventually 
applied to three-dimensional wing code theoretical thrust 
coefficients for a surface that also has zero thickness. 

Shown in figure 4 is an example of the variation of 
the attainable and theoretical thrusts with angle of attack 
for a given normal airfoil section at a given normal Mach 
number. Inset sketches show pressure distributions for 
6°, 12°, 18°, and 24° angles of attack. As figure 4 shows, 
the limitation imposed by vacuum pressures can be quite 
severe for high angles. In figure 4, the attainable thrust 
factor K, is shown as a function of angle of attack; how- 
ever, as discussed previously, the connection between 
two-dimensional airfoil sections and the three- 
dimensional wing sections is made through theoretical 
leading-edge thrust coefficients and not angle of attack. 
To find an appropriate three-dimensional wing section 
angle of attack to match a two-dimensional section angle 
of attack would be difficult, if not impossible, because of 
the extreme variation of upwash just ahead of the wing 
leading edge. The theoretical thrust coefficients provide 
a better connection because of the dependence of these 
coefficients on linearized theory singularity strength, 
which is a measure of pressure levels in the vicinity of 
the section leading edges. When pressure limiting has 
only a small effect, as it does for low Mach numbers and 
low angles of attack, the subsonic airfoil computer code 
gave values of theoretical thrust greater than that for a 
zero-thickness airfoil. Thus, K t can be greater than 1.0 
with a maximum value that tends to increase with 
increasing airfoil thickness. Because experimental data 
show little or no evidence of the theoretical benefit of air- 
foil thickness on attainable thrust given by the two- 
dimensional airfoil computer code, the attainable thrust 
factor K t , as shown in figure 4, is restricted to values of 
1.0 or less. 2 


2 In retrospect, an alternative procedure could have been applied. 

An attainable thrust factor defined as the ratio between thrust coef- 
ficients with and without pressure limiting ((2jt sin 2 a)/p„ replaced 
by c* for M n = 0) would automatically limit K t to values less 
than 1 .0. Although this alternative procedure has some attractive 
features, the resultant method would not be expected to give signif- 
icantly different results. As shown later, experimental data are used 
to calibrate the method. A different calibration would compensate 
for changes in the K, factor. 


After the vacuum pressure-limited thrust coefficient 
data are determined for the wide range of airfoil sections 
depicted in figure 3, the next step is to represent the data 
by empirical equations for use in automated calculations. 
The representation process is quite involved and was 
developed after considerable trial and error. For the inter- 
ested reader, a discussion of the strategy employed is 
given in appendix A. As discussed in appendix A, the use 
of a theoretical thrust parameter P tt and an attainable 
thrust parameter ^at provided the means of incorporating 
a range of airfoil geometric properties in a simplified rep- 
resentation of the attainable thrust factor K t . Results of 
the data representation are shown in figure 5 in the form 
of P at given as a function of P ttm Each of the five plots in 
figure 5 shows results for a given value of the radius 
index. 

The curves shown in figure 5 represent a fairing 
of the data provided by a single equation derived in 
appendix A to cover Mach numbers ranging from 0 to 
nearly 1 .0, maximum thickness ratios from (t /c) n = 0 
to 0.15, locations of maximum thickness from T] = 0.1 
to 0.5, and leading-edge radius indices from 0 to 1.2. The 
equation is 



or 


K, 




f x l f 11 Y 1 

«2' 

C ( ,A 



0.09 



_ 



with K t limited to values less than 1.0 where 


( 1 ) 



e l 

e 2 

*3 


0- 4C-0.7 

1- C - 3 ' 0 

-0.32^° -0.3 


The limitation of K t to values no greater than 1.0 permits 
attainable thrust to equal, but not exceed, theoretical 
thrust values defined by lifting surface theory. In the 
curve-fitting exercise, primary attention was given to 
representation of factors near the middle of the P at range. 
In addition, the greatest emphasis was placed on data 
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representing nominal airfoil parameters of (x/c) n = 0.09 
and rj = 0.5. Thus, the system is less accurate for extreme 
airfoil shapes, particularly for thin airfoils with forward 
locations of maximum thickness and sharp or nearly 
sharp leading edges. 

In figure 5, the decrease in ^at with increasing P n is 
clearly shown, as is the strong dependence on Mach 
number. A comparison of the plots in figure 5 shows the 
effect of increasing leading-edge radius. A sharp leading 
edge (fig. 5(a)) produces a substantial level of attainable 
thrust. Through an oversight, the previous attainable 
thrust method of reference 1 did not account for any of 
this thrust. For a leading-edge radius of zero, the pre- 
vious method gave a thrust of zero. Sketch F is an exam- 
ple of the variation of attainable thrust with increasing 
leading-edge radius for a 9-percent thick airfoil with 
maximum thickness at the 50-percent chord station at an 
angle of attack of 12° at M n = 0.5. 



0 .3 .6 .9 1.2 


Sketch F 


In figure 5 the considerable dependence of the 
attainable thrust parameter on maximum thickness and 
its location is not clearly evident. Sketch G shows the 
variation of attainable thrust with increasing thickness 
for an airfoil with a leading-edge radius index of 0.3 and 
maximum thickness at the 50-percent chord station at the 
same flow conditions. The nearly linear dependence of 
thrust on thickness clearly illustrates the importance of 
thickness and its frontal projected area in the develop- 
ment of thrust. 

The effect of the location of maximum thickness on 
developed thrust is illustrated in sketch H for the same 
nominal conditions of maximum (x/c)„ = 0.09, r i n = 0.3, 
a = 12°, and M n = 0.5. Benefits of a more forward loca- 
tion of projected frontal areas on which thrust is devel- 
oped are clearly shown. However, these thrust benefits 
are achieved at the expense of a tendency toward 
increased profile drag for such sections. 



0 .05 .10 .15 


(Vc) n 
Sketch G 



0 .2 .4 .6 


Sketch H 

Equivalent Mach Number Concept 

fr '/j.Iim £p,vac) 

Equation (1) was developed to account for the reduc- 
tion in attainable leading-edge thrust resulting from the 
application of realistic constraints on local pressure coef- 
ficients. A limiting pressure defined by the vacuum pres- 
sure coefficient has been shown to have a powerful effect 
on the amount of theoretical leading-edge thrust that can 
actually be realized. However, even more severe limits 
on achievable thrust are experienced in the real flow over 
airfoil sections when the local flow lacks sufficient 
energy to negotiate turns about the airfoil surface without 
becoming detached from that surface. Establishment of 
values for these more severe limitations is addressed in 
the following section of this paper. Before that, a means 
of application of equation (1) to the estimation of attain- 
able thrust for values of limiting pressure other than the 
vacuum limit is developed. Equation (1) can be used for 
a full range of limiting pressures between 0 and C p vac 
by substitution of a properly defined equivalent Mach 
number M e for the normal Mach number M n . The substi- 
tute Mach number is defined by the following logic. As 
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Sketch I 


illustrated by the pressure distributions shown in sketch I 
for a given airfoil section at the same positive angle of 
attack, the pressure coefficient at any point on the airfoil 
will vary with Mach number according to the Prandtl- 
Glauert rule. Thus, if the limiting pressure C pjlim also 
changes in accordance with the Prandtl-Glauert rule 

Jl-M 2 

f~K 


the attainable thrust factor K t will be the same at all 
Mach numbers because both c* n and c f will have the 
same Mach number dependence. Then with M e selected 
so that lim (M e ) = Cp Vac (M e ), the appropriate value of 
K t for the normal Mach number under consideration 
is calculated by substitution of M e for M n in equa- 
tion (1). The required M e is determined by setting 
C p vac (M e ) = C p i im (M e ), the intersection point of the 
curves shown in sketch I, and by solving for M e . Thus, 


Qvac W e ) = C,(M J 


2 

,.2 

yM e 






2 

n 


f 



and after the solution for the equivalent normal Mach 
number. 



Experimental Two-Dimensional Airfoil Analysis 

(Cp iinj Calibration) 

To define practical values of the limiting pressure 
coefficient, an incomplete version of the present method 
was applied to experimental two-dimensional airfoil data 
for symmetrical sections. (See refs. 10-18.) Correlations 
of axial-force coefficients predicted by this incomplete 
present method with experimentally determined axial- 
force coefficients, as shown in the examples of figure 6, 
were used to determine, by trial, values of C p j im that 
would match the experimental trends. For these symmet- 
rical sections, A c A is simply the negative of c* n . The 
example correlations in figure 6 were chosen to represent 
the procedures that were applied to the large amount of 
data available in the references. These data had a range 
of airfoil maximum thicknesses from 4 to 15 percent of 
the chord and locations of maximum thickness from 10 
to 42 percent of the chord. However, leading-edge radius 
indices had only a small range of 0.24 to 0.33. Mach 
numbers ranged from 0.03 to 0.90, and Reynolds num- 
bers varied from less than 1 x 10 6 up to 30 x 10 6 . 

To evaluate limiting pressure coefficients, equa- 
tions (1) and (2) were combined in a computational pro- 
cess in which K t and c* n were calculated for a series of 
trial Cp jin, values to find the value that most closely 
matched the experimental data. In matching the trial 
curve fit to the experimental data, particular attention is 
given to breakaway of the experimental axial force from 
the theoretical leading-edge full-thrust curve. For most of 
the plots, this breakaway point can be established with 
reasonable certainty. For other plots a breakaway point 
of the experimental data is not readily obvious. The prob- 
lem occurs because axial-force data were not presented 
directly for some of the experimental investigations; 
the axial force had to be derived from lift- and drag- 
coefficient data. For some of these data, as in the exam- 
ples in figure 6(c), C D data were omitted even though C L 
versus a data were shown in the plots. Generally, this 
occurred for angles of attack in excess of the angle at 
which maximum lift was obtained. Presumably the 
authors wanted to avoid two values of drag at a given lift 
coefficient. Because of the absence of Q> data, axial 
force could not be calculated for all of the angles of 
attack listed in the reference data. However, extrapola- 
tion of the C D curve subject to nearly maximum C L or 
subject to c D values greater than the maximum plot 
values would lead to C A values on or above the curves 
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corresponding to the specified limiting pressure coeffi- 
cients. Thus, when properly interpreted, these data also 
provide estimates of limiting pressures. For the lower 
plots in figure 6(c), complete data were provided for the 
highest angle of attack shown. These experimental data 
closely follow the theoretical full-thrust curve for the 
entire plot and show no breakaway tendencies. For this 
case, only the fact that limiting pressures are less than 
(i.e., greater negative value) those that would cause a 
breakaway at the highest angle of attack at which data are 
shown could be determined. The examples shown in fig- 
ure 6 include curves of calculated A C A for the selected 
value of C p >lim and for values of ±20 percent of C p lim . 
Obviously, establishment of precise values of limiting 
pressure coefficients is not possible; however, as shown 
in figure 6, attainable thrust and, especially, its break- 
away point can be predicted with reasonable accuracy 
with only approximate limiting pressure estimates. 

Data from figure 6 and from a large number of other 
plots of axial-force coefficient versus angle of attack 
(87 plots in all) were used in figure 7 to define the depen- 
dence of the limiting pressure coefficient on Mach num- 
ber and Reynolds number. Although the C p j im data 
showed some dependence on airfoil section geometric 
characteristics, this dependence was not systematic and 
was small compared to that of M n and R n . The limiting 
pressure is shown as a function of Mach number for a 
series of nominal Reynolds numbers. Nominal Reynolds 
number data were obtained by either interpolation or 
extrapolation of Reynolds number data such as that 
shown in sketch J. All of the data shown in figure 7 were 
obtained for relatively smooth airfoil sections with no 
attempt to establish turbulent flow through the use of 
artificial roughness. Thus, the Reynolds number effect 
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shown here includes the effect of transition from laminar 
to turbulent flow. A large amount of scatter is noted in 
the data, especially for the nominal Reynolds number of 
3.0 x 10 6 . At this Reynolds number, evidence of transi- 
tional flow can reasonably be expected. The solid curves 
shown in figure 7 represent an attempt to mathematically 
express the variation of the limiting pressure coefficient 
with Mach number and Reynolds number. The curve fit 
was chosen to closely represent the upper limit of the 
data. Thus, a system that uses this empirical relationship 
to predict attainable thrust tends to represent idealized 
conditions, i.e., smooth surfaces and stable uniform flow 
with little turbulence ahead of the airfoil or wing. 

The curve-fit equation chosen to represent the data 
has the form 


"p, lim 


Y M n 


f Rx 10 6 


\R n x 10 + K j 


( 3 ) 


which was chosen to give a limiting pressure approach- 
ing the vacuum limit for infinitely large Reynolds num- 
bers and a limiting pressure approaching 0 for infinitely 
small Reynolds numbers. Values of K and the exponent 
<? 4 were chosen by a trial-and-error process to closely 
match an upper bound of the empirically derived C p lim 
distribution. These values are 


K = 


10 


8(1 ~M„) 


-0 75 

e 4 = 0.028M n 

The data and the curve fits in figure 7 show a strong 
dependence of the limiting pressure on Mach number 
equal, at least, in importance to the dependence on 
Reynolds number. The tendency for the limiting pressure 
to decrease with increasing Mach number was expected 
because C p vac decreases with increasing M n . The exist- 
ence of a limiting pressure peak and the drop-off for 
Mach numbers less than about 0.15 were not anticipated. 
To search for validation of this Mach number trend and 
to further explore Reynolds number effects, the Navier- 
Stokes computer code described in reference 19 was 
employed as described in appendix B. Basically, the 
computer code was used to create two-dimensional air- 
foil axial-force data, which was treated in the same fash- 
ion as the experimental data, to define another set of 
limiting pressures. Because of the time-consuming 
nature of the computer code, only a few representative 
cases were handled. The Navier-Stokes computer code 
was first used to provide data for a NACA 0009 airfoil at 
a Reynolds number of 9 x 10^ for a series of Mach num- 
bers. For these conditions shown in figure 7(d), C p y m , as 
defined by the Navier-Stokes computer code data, has a 
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variation with Mach number si mi lar to that given by the 
data but at a slightly higher level. This is understandable 
because the airfoil model can be expected to have surface 
roughness not modeled in the theoretical airfoil shape, 
and the wind tunnel flow could have turbulence levels 
higher than anticipated. Even with this additional confir- 
mation, some uncertainty remains about the evaluation of 
limiting pressure for normal Mach numbers below 0.1. 
Thus, results of the computer codes employing this 
attainable thrust estimation technique are questionable 
for Mach number and sweep angle combinations in 
which M cos Aj e <0.1. For example, results for a 70° 
swept leading-edge wing could be somewhat suspect for 
ffee-stream Mach numbers less than about 0.3. 

Navier-Stokes computer code data depicting the 
dependence of limiting pressures on Reynolds number 
are compared with results given by the curve fit 
in sketch K. The Navier-Stokes results are for an 
NACA 0009 airfoil at a Mach number of 0.16. For this 
Mach number, the Navier-Stokes results indicate a more 
beneficial effect of increasing Reynolds number than is 
given by equation (3) used in the present method. Fur- 
ther, more complete studies using Navier-Stokes com- 
puter codes in conjunction with appropriate experimental 
verification may eventually lead to a more accurate rep- 
resentation of limiting pressures. In particular, a need 
exists for better definition of limiting pressure coeffi- 
cients at high Reynolds numbers and low Mach numbers. 
The use of such additional data in a revised calibration of 
the present attainable thrust prediction method is dis- 
cussed in the section entitled “Future Limiting Pressure 
Redefinition With Use of AER02S Computer Code.” 
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omitted. The purpose is to show in a single plot the influ- 
ence of both Mach number and Reynolds number on lim- 
iting pressure trends. At high Mach numbers, the limiting 
pressure obviously is relatively insensitive to Reynolds 
number and is severely constrained by the vacuum pres- 
sure limit. However, at low Mach numbers, a strong sen- 
sitivity to Reynolds number exists. With the present 
curve fit, the greatest changes with Reynolds number 
occur for a Mach number of about 0.1. However, this 
strong sensitivity does not lead to correspondingly large 
changes in leading-edge thrust because C p ;lim acts only 
as a bound in the pressure integral for c* n . 

These limiting pressure data from the analysis of 
two-dimensional airfoil experimental data (eq. (3)) are 
used to define an equivalent Mach number M e (eq. (2)) 
for use in the theoretical two-dimensional airfoil analysis 
(eq. (1)) to estimate the dependence of attainable thrust 
for two-dimensional airfoils on airfoil geometric proper- 
ties and flow conditions. In essence, the experimental 
two-dimensional airfoil data provide the calibration of 
the system. The product of this two-dimensional analysis 
is a factor K, defined by equation (1) that relates the 
attainable two-dimensional thrust to the theoretical two- 
dimensional thrust by 


Guidelines for Present Method Application to 
Three-Dimensional Wing Computer Codes 

The present method for the prediction of attainable 
leading-edge thrust described in this paper is intended to 
be incorporated into linearized theory lifting surface 
computer codes (e.g., AER02S and WINGDES), which 
provide a calculation of theoretical leading-edge thrust. 
The first step in programming of the attainable leading- 
edge thrust method for inclusion in linearized theory lift- 
ing surface computer codes is implementation of the sim- 
ple sweep analysis to define the equivalent normal or 
two-dimensional airfoil geometric characteristics and 
flow conditions for a series of streamwise sections from 
the wing root to the wingtip. For each of a large number 
of wingspan stations, the following terms define the nor- 
mal section geometric characteristics: 


Sketch K 


C n = c cos A le 


The variation of the limiting pressure coefficient 
with Reynolds number is depicted in figure 8. The curve- 
fit data shown here were obtained from equation (3). The 
experimental data used in generation of the curve fit are 


'*) = 
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cos A le 


which accounts for differences between the more realis- 
tic limiting pressure coefficients and the vacuum pres- 
sure coefficient. 



In addition, for each span station, the normal Mach 
number, the normal Reynolds number, and the normal 
section thrust coefficient are defined as follows: 


The attainable thrust factor K t is calculated from 



with K t limited to values less than 1.0 where 


M„ = M cos A le 


R n = R J cos A le 
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or with c t nondimensionalized as in AER02S and 
WINGDES 
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The theoretical thrust coefficient c t is supplied by the 
linearized theory lifting surface analysis. This coeffi- 
cient accounts for variations of the theoretical leading- 
edge full-thrust coefficient with such factors as free- 
stream Mach number, wing planform, and wing twist and 
camber. 

For each span station, the limiting pressure coeffi- 
cient is then calculated from 


"p, lim 
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where 
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The limiting pressure coefficient is then used to calculate 
the equivalent Mach number 
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The ratio of attainable theoretical thrust K t for each 
normal airfoil section is then applied to the theoretical 
thrust for the wing at the wing spanwise station from 
which the normal section was derived to arrive at a span- 
wise distribution of attainable thrust. The attainable 
thrust section coefficient is then 


where c t is the theoretical leading-edge full-thrust coeffi- 
cient defined by the attached flow linearized theory lift- 
ing surface solution. 

For a given wing and a given set of test or flight con- 
ditions, the theoretical leading-edge thrust c t given by a 
linearized theory computer code will vary with span 
position and wing angle of attack. The attainable thrust 
factor K t determined by the attainable thrust method will 
vary with those quantities and with other factors includ- 
ing local leading-edge sweep angle and wing section 
characteristics. The effect of wing twist and camber or 
flap deflection is accounted for in the calculated c t value. 
Because the attainable thrust factor is dependent on the 
theoretical thrust, it is also influenced by wing twist, 
wing camber, or flap deflection. The same spanwise inte- 
gration techniques used in the lifting surface computer 
code can be employed to calculate wing attainable thrust 
coefficients. 
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Present Method Validation 

To illustrate the applicability of the present attain- 
able thrust prediction method, a series of comparisons 
of theoretical and experimental data for wing-body 
configurations is presented in figures 9-18. The 
AER02S computer code described in reference 7 and the 
WINGDES computer code described in reference 8 were 
modified to replace the previous attainable thrust predic- 
tion method with the present version. For some exam- 
ples, results using both versions are shown. Computer 
code-estimated forces include attainable thrust and the 
effects of a separated vortex whose strength is deter- 
mined by the Polhamus leading-edge suction analogy 
and whose location is given by delta wing empirical data 
(computer code vortex option (1)). Computer code 
results are designated by the solid line for the present 
method and by the long-dash-short-dash line for the pre- 
vious method. Computer code results without the addi- 
tion of either attainable leading-edge thrust or vortex 
forces are shown by the short-dash line. Because some 
amount of either of these forces is almost always present, 
this curve gives a conservative estimate of the drag upper 
limit. As an additional reference, a theoretical leading- 
edge full-thrust curve given by the computer code is 
shown on the axial-force coefficient plot and on the lift- 
drag plot. The purpose of the attainable thrust prediction 
is to provide an estimate of the portion of the theoretical 
thrust that can actually be achieved. These reference 
curves help in an assessment of the success of the 
prediction. 

Because the purpose of this attainable thrust study is 
the development and assessment of methods for predic- 
tion of drag associated with the generation of lift, the 
zero-lift drag coefficient 0 for a flat wing configura- 
tion (a wing with no twist or camber) used to construct 
the theoretical curves was obtained from experimental 
data. This was accomplished by setting Q>, o equal to 
the limiting value of the axial-force coefficient for a 
flat wing configuration as a approaches 0 (C D = C A at 
a = 0°). Use of the C A rather than the C D curves gives a 
more accurate result because of the lesser sensitivity to 
angle of attack or C L variation. 

Figure 9 shows data at a Mach number of 0.4 and a 
Reynolds number of 2.5 x 10 6 for a 60° swept delta wing 
fighter with a rounded leading-edge NACA 64A00X air- 
foil. Correlations using this experimental data from refer- 
ence 20 were originally shown in reference 6. For the 
undeflected flap or flat wing case in figure 9(a) and the 
20° deflection case in figure 9(b), the present attainable 
thrust method predicts thrust that is slightly greater (i.e., 
more negative C A ) than the previous method but nearly 
equal to the experimental values. 


For these data and for all of the theoretical- 
experimental correlations to be shown, little or no differ- 
ence is noted in predicted normal-force coefficients and 
pitching-moment coefficients because of the change in 
the attainable thrust prediction. Differences between 
experimental and theoretical C N and C m curves shown in 
figure 9 are not addressed here. However, the C N differ- 
ences are responsible for the lack of improvement in C D 
prediction for the present method in spite of the C A 
improvement. Note that for the deflected flap case, the 
difference between no-thrust and full-thrust theoretical 
C A curves is smaller at a given positive angle of attack. 
Thus, the deflected flap configuration is less sensitive to 
the attainable thrust prediction than is the undeflected 
flap or flat wing configuration. 

Data for the same 60° swept delta wing fighter at a 
higher Mach number of 0.8 are shown in figure 10. Here, 
the newer attainable thrust prediction is slightly poorer at 
angles of attack above 10°. The prediction of drag does 
not change below lift coefficients of about 0.5 and 
changes by only a small amount up to C L = 0.8. This 
poorer prediction of wing characteristics for high-lift 
coefficients as sonic speeds are approached is to be 
expected. The AER02S and WINGDES computer codes 
are not applicable at transonic speeds. 

Figure 1 1 shows data at a Mach number of 0.4 for a 
wing-body configuration identical in all respects to that 
in figure 9 except for the substitution of a sharp leading- 
edge airfoil. The leading edge of this airfoil, which can 
not practically have a sharp leading edge, was estimated 
to have a leading-edge radius of 0.002 in. With this 
radius, the present method gives a very good estimate of 
the thrust actually achieved. Sketch L illustrates an 
important point concerning the sensitivity of thrust to the 
leading-edge radius in the previous and present system. 
With the previous method, a wide variation between 
results with zero radius and the small estimated radius 
was obtained. For the present method, that range is 
reduced, and even if the radius were zero, a significant 
amount of thrust would remain. Thus, the present method 
is not nearly so dependent on the skill of the user in 
choosing an appropriate radius for a particular case. 

Data for a 44° swept trapezoidal wing fighter at a 
Mach number of 0.4 and a Reynolds number of 1.9 x 10 6 
are shown in figure 12. These data from reference 20 
were also used in correlations presented in reference 6. 
Little difference between the predictions given by the 
two methods is seen up to an angle of attack of about 8° 
or a C N of about 0.5, and both methods agree well with 
the experimental data. However, beyond this angle of 
attack, there is evidence of a severe flow breakdown 
which the theory cannot handle. The reduction in 
leading-edge sweepback angle from that of the previous 
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Present method 



Sketch L 


correlations may be responsible, at least in part, for the 
earlier and more severe breakaway of the experimental 
axial-force coefficient from the computer code predic- 
tion. For large sweep angles, drastic flow separation 
would be expected to begin near the wingtip at relatively 
small angles and progress inward as the angle of attack is 
increased. This would result in a thrust loss that is more 
gradual than that for small sweep angles where drastic 
flow separation would be expected to occur more 
abruptly. As noted earlier, axial-force breakaway can be 
very severe for two-dimensional airfoils (zero sweep 
angle) at low Mach numbers. 

When leading- and trailing-edge flaps for the trape- 
zoidal wing fighter are deflected to 20°, the results 
shown in figure 12(b) are obtained. Again, the present 
method predicts a somewhat greater amount of attainable 
thrust but fails to adequately predict the experimen- 
tal values of achieved thrust. At a lift coefficient of 
about 0.6, where the deflected flaps provide a substantial 
drag reduction, both the present and the previous meth- 
ods give a reasonable prediction of the performance. 
Again, the change in attainable thrust prediction with 
deflected flaps has less effect than that with undeflected 
flaps. Because of this smaller effect for wings with 
deflected flaps (or twisted and cambered wings), the 
remainder of the theoretical-experimental correlations 
will be restricted to flat or slightly cambered wings. 

Data for a supersonic transport configuration from 
reference 21 are shown in figure 13. Comparisons of this 


data with the previous method were given in reference 6. 
The present method gives a better prediction of C A and 
C D . Experimental axial-force data for this wing with very 
high leading-edge sweep angles do not display the sharp 
upturn away from the computer code estimate shown in 
the previous example. For this wing at the specified test 
Mach number of 0.25, the normal Mach number is about 
0.08 for the inboard panel and about 0.12 for the out- 
board panel. The reasonably good correlation of com- 
puter code results with the experimental data should help 
to allay some previously expressed concerns regarding 
the proper selection of limiting pressures for very low 
normal Mach numbers. Similar results, which are not 
included in this collection, were obtained for another 
supersonic transport configuration tested in a different 
facility. 

Data obtained at a supersonic Mach number for a 
wing with three different airfoil sections are shown in 
figure 14. These data from reference 22 were used in cor- 
relations presented in reference 8. For all leading edges, 
which included the sharp leading edge (fig. 14(a)), the 
standard leading edge (fig. 14(b)), and the blunt leading 
edge (fig. 14(c)), little discernible difference is noted 
between the predictions given by the two methods. Both 
give results in good agreement with the experimental 
data. For the sharp leading edge, the previous method 
with an estimated leading-edge radius works surprisingly 
well. However, as noted previously, it is much more sen- 
sitive to the radius selection than the present method. 
Without the assumed radius, the previous method would 
give a C A value of 0.01 independent of the angle of 
attack. 

An experimental program designed to examine the 
effects of wing leading-edge radius and Reynolds num- 
ber on wing performance was reported in reference 23. 
Some examples from that report, which show the effects 
of leading-edge radius, are shown in figure 15. For these 
comparisons, only the present version of the attainable 
thrust method was used. The configuration that was 
selected had a leading-edge sweep angle of 61.7° and an 
unswept trading edge. The test Mach number was esti- 
mated to be 0.22. The Reynolds number for the data pre- 
sented in figure 15 has a small range between 9.72 x 10 6 
and 10.14 x 10 6 . For the sharp leading edge (fig. 15(a)) 
a very small leading-edge radius of 0.002 in. was 
assumed as before. For the other sections, leading-edge 
radii specified in the report were used. In the three plots 
in figure 15, the leading-edge radius has a significant 
effect on performance, and the present attainable thrust 
method predicts this performance quite well. The only 
appreciable discrepancy between theory and experiment 
occurs for the sharp leading-edge airfoil section. As 
noted in discussions of the present method development, 
the method is expected to be least accurate for thin wings 
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with sharp leading edges and forward locations of 
maximum thickness. The wing section thickness ratio is 
only 0.02 at the wing-body juncture, and the location of 
maximum thickness ri is as low as 0.12 for much of the 
wingspan. 

Figure 16 shows data for the same configuration as 
in figure 15 with the small leading-edge radius (actually 
more accurately characterized as a standard leading-edge 
radius with r i ~ 0.27) at three additional Reynolds num- 
bers. At the lowest Reynolds number (fig. 16(a)), the 
attainable thrust is somewhat overestimated, and at the 
highest Reynolds number (fig. 16(c)), it is slightly under- 
estimated. However, in general, both the axial force 
and the drag are predicted reasonably well at all three 
Reynolds number. 

The data presented in figures 17 and 18 provide 
additional verification of the present attainable thrust 
prediction method and also illustrate a problem that is 
sometimes encountered in assessing wing drag-due-to- 
lift characteristics from wind tunnel data. Data for a 44° 
sweptback wing of aspect ratio 2 with a 3-percent-thick 
rounded nose section in Mach number 0.61 flow taken 
from reference 24 are compared with theoretical results 
in figure 17. In figure 17(a), for a Reynolds number of 
1.9 x 10 6 , rather peculiar behavior of the experimental 
C A occurs in the angle-of-attack range from ±3°. This 
test model, which employed no artificial means of stabi- 
lizing the boundary layer, may have generated a flow 
with large laminar flow regions at and near 0° angle of 
attack. With increasing angle of attack, the flow appar- 
ently became predominantly turbulent. The assessment 
of drag-due-to-lift characteristics requires the establish- 
ment of a Or 0 value reflecting the type of flow (laminar- 
turbulent transition) present at lifting conditions. The 
experimental C D (and C A ) at a = 0° of about 0.0070 does 
not represent such a value. However, as shown in the C A 
plot, an appropriate C D 0 can be found by fitting the 
curve given by the computer code to experimental axial- 
force coefficients in the range where the slopes agree. 
The C O 0 value of 0.0100 given by this process, when 
combined with the computer code drag-due-to-lift pre- 
diction, provides good correlation with the measured 
results up to about 12° angle of attack. However, a sig- 
nificant discrepancy in the pitching moment, which is 
more than that found for most of the other configura- 
tions, is shown. 

Data for the 44° sweptback wing at a higher 
Reynolds number of 4.8 x 10° are shown in figure 17(b). 
Here, the application of the curve-fit technique also 
requires the establishment of a representative Q>, 0 value. 
However, in this case, only a single C A data point (that 
closest to a = 0°) departs significantly from the theory 
curve. Again, the drag polar data are well predicted. 


Figure 18 presents data for the same configuration as 
in figure 17 with the exception of the substitution of a 
sharp leading edge for the rounded leading edge. For the 
calculations a constant leading-edge radius of 0.002 in. 
was assumed. Only a small Q>, 0 correction was required 
for these test data. Comparison of these results with those 
in figure 17 shows almost no penalty in either theoretical 
or experimental C A and c D associated with the sharper 
leading edge. A small loss in leading-edge thrust is offset 
by a lower Q>, 0 - 

Sketch M helps in a discussion of the importance of 
a proper establishment of C D 0 in assessing drag-due-to- 
lift characteristics. The data shown here are for the con- 
figuration analyzed in figure 17. The circles represent the 
experimental data, and the solid lines a possible fairing 
of that data. The drag at zero lift displays an increase 
with increasing Reynolds number rather than a decrease 
as would be expected. This is in marked contrast to the 
behavior of the drag of this same configuration at higher 
lift coefficients (e.g., C L = 0.4). The drag at zero lift is 
likely influenced by flow changes (transition from lami- 
nar to turbulent) not experienced at higher lift coeffi- 
cients. At higher Reynolds numbers where the boundary 
layer at both lift coefficients is likely to be predominately 
turbulent, the two curves would tend to be parallel. An 
assessment of the required C D ,o correction for a given 
wing, as represented by the arrow in sketch M, can be 
determined by use of the C A curve-fitting technique pre- 
viously discussed. 
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Sketch M 


A review of all the correlations presented in fig- 
ures 9-18 shows that, in general, the initial breakaway of 
axial force from the theoretical full-thrust curve is pre- 
dicted quite well. However, for a wing with a 44° 
leading-edge sweep angle at a Mach number of 0.4, a 
more severe loss of thrust occurred beyond this point 
than that predicted by the present attainable thrust 
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method. For high-sweep angles (A le = 60° or more), the 
axial force is predicted with reasonable accuracy for 
angles of attack up to about 16° or 20°. Highly swept 
wings are expected to display a more gradual onset of the 
effects of drastic flow separation because the initial sepa- 
ration tends to occur near the wingtip and progress 
inward as the angle of attack is increased. This is in con- 
trast to more moderately swept or unswept wing sections 
where the drastic separation can occur all at once. 

Future Limiting Pressure Redefinition With 
Use of AER02S Computer Code 

As explained in the section of this paper enti- 
tled “Experimental Two-Di m ensional Airfoil Analysis 

jjjjj Calibration),” the effects of Reynolds number and 
Mach number on attainable thrust are derived empirically 
by use of available two-dimensional airfoil experimental 
data. Thus, the accuracy of the whole system depends on 
the availability and accuracy of these experimental data. 
A new set of two-dimensional airfoil experiments tai- 
lored specifically for this problem would make a signifi- 
cant contribution to the technology. In addition, as 
indicated in a preliminary study included in the method 
development section, two-dimensional Navier-Stokes 
solutions may provide valuable additional information 
including extension of the capabilities of the experimen- 
tal approach. In particular, better definition of limiting 
pressure coefficients is required at high Reynolds num- 
bers for low Mach numbers. In the hope of encouraging 
such further activities, some suggestions for a revised 
calibration of the system are given. 

The problem in calibrating the attainable thrust 
method is to find, by iteration, C p;lim values that give 
axial-force coefficients fitting the experimental data as it 
breaks away from the theoretical full-thrust curve. To 
facilitate any future recalibration of the system, the 
AER02S computer code has been modified to provide a 
two-dimensional airfoil solution. To use the computer 
code in the two-dimensional mode, set JBYMAX = 1 
and input data for two identical wing sections at the root 
and tip stations of a rectangular wing of any span. In 
this mode the computer code can be used to supply addi- 
tional calibration data. This is accomplished by input of 
xmcplt values other than 1.0, which act as multipliers 
of the limiting pressure as now defined. When a chosen 
xmcplt value gives a A C A curve matching the experi- 
mental data for a range of angles of attack near the 
breakaway point, the corresponding limiting pressure 
coefficient included in the output data is taken to be the 
value for the input Mach number and Reynolds number. 
A collection of such data for a large number of airfoils 
and test conditions then can be used to generate a revised 


curve-fit equation that might, but not necessarily, be sim- 
ilar in form to that presently used. 

The part of the AER02S and WINGDES computer 
codes that is to be changed to accommodate a new for- 
mulation of the limiting pressure coefficient definition 
can be found by conducting an automated search for 
_cplt_. Inclusion of the space before and after the 
name excludes other extraneous lines. In the AER02S 
computer code, the _cplt_ definition occurs twice, 
once for each surface. 

All of the two-dimensional experimental data 
employed in the present analysis were for symmetrical 
sections. Airfoils with camber can also be used if neces- 
sary. For such airfoils, ordinate tables for a mean camber 
surface (halfway between lower and upper surface ordi- 
nates at a given chord station) are required. The com- 
puter code then provides an axial-force variation with 
angle of attack differing from the symmetrical section 
values, but the iteration process is no different. The goal 
is still to find C p lim values by iteration that match com- 
puter code and experimental axial-force data as the 
experimental data breaks away from the theoretical 
leading-edge full-thrust curve. 

Concluding Remarks 

A revised method (called the present method herein) 
for predicting attainable leading-edge thrust and its effect 
on airplane aerodynamic performance has been devel- 
oped. The present method accommodates a greater range 
of airfoil shapes from very sharp to very blunt leading 
edges. It is also based on a wider range of Reynolds num- 
bers than the range of the previous method. An additional 
aesthetic appeal of the present approach is a relationship 
between flow over the wing surface and the correspond- 
ing two-dimensional airfoil analysis that is more consis- 
tent with simple sweep theory concepts and accounts 
more accurately for the actual flow over the wing. The 
present method, when employed in the aerodynamic 
design and analysis computer codes WINGDES and 
AER02S, provides results that generally give improved 
correlation with experimental wing-body axial-force data 
and reasonable estimates of the measured drag. Applica- 
tion of the present method for sharp leading-edge airfoil 
sections is simplified because the uncertainty concerning 
selection of an appropriate leading-edge radius to replace 
the nominal value of zero no longer exists. The present 
method can be expected to be most accurate for wings 
with standard airfoil sections and least accurate for non- 
standard sections, particularly thin sections with sharp or 
nearly sharp leading edges and forward locations of max- 
imum thickness. 
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Provisions have been made for a simplified recali- 
bration of the empirical part of the present method when 
experimental two-dimensional airfoil data for a wider 
range of airfoil shapes or test conditions become avail- 
able. Particularly, two-dimensional data are needed at 
high Reynolds numbers and low Mach numbers for low- 


speed flight conditions of large highly swept wings such 
as those employed on supersonic transport aircraft. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
November 30, 1995 
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Appendix A 

Strategy for Thrust Factor Data 
Representation 

The process used to collapse the large amount of 
two-dimensional attainable thrust factor data generated 
by use of the airfoil computer code to the relatively nar- 
row bands shown in figure 5 is described. In addition, the 
development of an empirical equation to represent these 
bands and the dependence of the attainable thrust factor 
on the Mach number and airfoil geometric characteristic 
terms used in figure 5 is discussed. 

For a given leading-edge radius index and Mach 
number, the collapse of data for the wide range of airfoil 
geometries shown in figure 3 is accomplished through 
use of the theoretical thrust parameter 



and the attainable thrust parameter 



Values of the exponents e t and e 2 are different for each 
radius index. The nominal value for the maximum thick- 
ness ratio is 0.09, and the nominal value for the location 
of maximum thickness is 0.5. The selected parameters 
collapse data for other thicknesses and locations of maxi- 
mum thickness on the data band defined by the nominal 
values. 

The attainable thrust parameter contains a section 
thickness term because the integrated attainable thrust 
coefficients given by the two-dimensional airfoil com- 
puter code show an almost linear increase with increas- 
ing thickness. Attainable thrust coefficients with no 
limiting (obtained from the M n = 0 data, C p ilim = — °°) 
give the maximum values. In sketch N, the maximum 
thrust factor K t is shown as a function of the section 
thickness. The data in the sketch cover radius indices 
r i n = 0, 0.15, 0.30, 0.60, and 1.20, as well as locations 
of maximum thickness r| of 0.1, 0.3, and 0.5. Most 
of the data cluster near a tine given as maximum 
K t = 1 + (x/c) n . This line is defined primarily by the 
data for a standard leading-edge radius r i n = 0.3 and for 
a nominal value of the location of maximum thickness, 
r| = 0.5. The data showing the greatest departure from 
this line are for the sharp leading edge r in = 0 and the 
most forward location of maximum thickness T| = 0. 1 . 
Because of the discrepancy shown here, and others noted 



0 ^ .05 .10 .15 


M n 

Sketch N 

elsewhere in this paper, the attainable thrust prediction 
method is least accurate for airfoils with sharp leading 
edges and forward locations of maximum thickness. The 
attainable thrust parameter definition used here aids in 
the collapse of data for the range of section thickness 
from 0.03 to 0.15. 

A sample of the strategy used to accomplish the 
reduction of the data to narrow bands is explained with 
the aid of sketches O and P. Here the object was to super- 
impose data for other thickness ratios on data for the 
nominal value of 0.09. Sketch O shows two-dimensional 
airfoil computer code results at a Mach number of 0.5 for 
an airfoil with a leading-edge radius index of 0.3 (stan- 
dard) and a maximum thickness location of 0.5. The 
attainable thrust parameter is shown as a function of the 
theoretical thrust c tr fi n . To simplify the sketch, data are 
shown for only three thickness ratios. 



Sketch O 
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Sketch P 


The purpose of the exercise is to find an exponent e 2 
for the theoretical thrust parameter that superimposes the 
data for thicknesses other than 0.09 on the 0.09 curve. 
The form of the transformation is taken to be 



so that 



e 2 = 


lo g( C f,»P„) 


1 = 0.09 

Jn 


lo g( c t,»P») 


31 


log 




Sketch Q 


nation resulted in the data organization shown in fig- 
ure 5. For a given radius index and Mach number, data 
for a range of thickness and a range of maximum thick- 
ness location now have only a small amount of scatter 
about values for the nominal thickness (x/c) n = 0.09 and 
the nominal maximum thickness location r| = 0.5. An 
ideal solution for this collapsing process would also pro- 
vide for superposition of the data for other Mach num- 
bers on the M n = 0.5 data in each of the plots in figure 5 
and, in addition, would allow all data for the different 
radius indices to be superimposed on the r i n = 0.3 data. 
Such a solution could not be found. Instead, a means of 
devising a single equation to provide a curve fit of the 
data in figure 5 was employed. 


The exponent e 2 is found from an examination of cross- 
plotted data such as that displayed in sketch P. Emphasis 
is placed on attainable thrust parameters in the middle 
of the range. For the example shown in sketch O, 
a thrust factor of 0.43 was chosen. For (x/c) n = 0.03, 
the corresponding c f „P n -is approximately 0.048, which 
requires an exponent of -1.57 to provide superposition. 
For (x/c) n = 0.15, the corresponding c tr fi n is approxi- 
mately 0.60, which requires an exponent of -1.55. The 
selected exponent is taken as an average from a number 
of such evaluations for other thickness ratios and a range 
of thrust factors. When the chosen value of e 2 = -1.60 is 
used, the separate curves of sketch O are reduced to 
essentially the single curve shown in sketch Q. 

A similar process was used to collapse the data for 
the other locations of maximum thickness onto the 
ri = 0.5 curves. The collapsing process for maximum 
thickness and location of maximum thickness in combi- 


The problem of relating the attainable thrust parame- 
ter to changes in Mach number and leading-edge radius 
as shown in figure 5 was handled by application of an 
equation of the form 

*Vt = ^t/ 3 

to the data bands through selection of the constant k and 
the exponent e 3 . When two points on the desired curve 
are specified, the exponent e 3 may be determined from 
the expression 

= l°g(P a t)l- lo g( P at) 2 
* 3 " log (F tt ) j - log (F tt ) 2 

The two points were chosen to represent the curve in the 
region of greatest interest, generally, near the ^at 
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midrange. For example, with r in = 0.3 (fig. 5(c)), the 
curve for M n = 0.5 is plotted to pass through the points 

[(P at )j = 0.60, (P tt ) 1 = 0.15] 

and 

[(P at ) 2 = 0.30, (P a ) 2 = 0.48] 


k = 0.14[1.0-(1.0- i yr. „)M„] + 0.iyr. n 


(1 -M n ) 0A8 < 1+ O 


to give an exponent e 3 ~ -0.60 and a constant k ~ 0.20. 
The exponent was found to have little dependence on the 
Mach number. Therefore, although k depends on both the 
Mach number M n and the radius index r in , e 3 depends 
only on the radius index. 

When the preceding process was carried out for each 
of the five Mach numbers for each of the five leading- 
edge radius indices, the following equations provided a 
reasonable representation of the constant and exponents: 




i * 010 a n 
e 2 = 1 ' 6r i,n ~ 3 - 0 


e 3 = -0.32r°;°-0.3 


The faired curves shown in figure 5 were defined by use 
of these constants and exponents in equation (1). 
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Appendix B 

Navier-Stokes Solutions for Two-Dimensional 
Airfoils 

An essential part of the attainable thrust prediction 
method is the determination of effective limiting pressure 
coefficients and their dependence on local Mach number 
and Reynolds number. As described previously, this is 
accomplished through the use of wind tunnel data for 
two-dimensional airfoil sections. However, recognize 
that wind tunnel test results are subject to some limita- 
tions. The results are dependent on the precision of the 
measuring devices, on model accuracy (including surface 
smoothness), and on wind tunnel flow quality, which can 
vary from test to test and from facility to facility. 
Furthermore, and perhaps more importantly, the informa- 
tion provided by currently available data does not cover a 
matrix of test conditions (Mach number and Reynolds 
number) as broad ranged and as closely spaced as 
desired. In particular, some uncertainty exists with regard 
to the dependence of limiting pressures on Mach num- 
ber for low Reynolds numbers. Also, a broader range of 
Reynolds numbers is required to provide a more confi- 
dent extrapolation to flight conditions. 

To address some of these concerns, an exploratory 
investigation of the use of a Navier-Stokes equation- 
solving computer code to provide a substitute for experi- 
mental data was performed. First, to validate the com- 
puter code for application to this problem, computer code 
solutions were compared with representative experimen- 
tal data. Then, computer code results were used to pro- 
vide some limited samples of the use of Navier-Stokes 
solutions to supplement or replace experimental data. 

The Navier-Stokes computer code CFL3D employed 
in this study is described in reference 19. The governing 
equations used in the two-dimensional version of this 
computer code are the two-dimensional, Reynolds- 
number-averaged Navier-Stokes equations written in 
conservative form and expressed in a body-fitted coordi- 
nate system. These equations are solved by an implicit, 
finite-volume, upwind spatially factored algorithm. More 
information about this computational method and its 
accuracy is contained in reference 19. Turbulent viscos- 
ity was simulated with the transport equation of Spalart- 
Allmaras. (See ref. 25.) 

A 2-D grid was constructed for each of the following 
airfoils: NACA 0006, NACA 0009, NACA 0012, NACA 
64-010 with 257 and 81 grid points in the stream wise and 
normal directions, respectively. Each grid was suffi- 
ciently clustered normal to the body surface to resolve 
the viscous boundary layer. The grids were also clustered 
at the airfoil nose and trailing edge. 


Computations were performed on the CRAY-2S at 
NASA Langley Research Center with the use of multi- 
block iteration to accelerate convergence. A converged 
solution for a typical case was obtained after approxi- 
mately 1500 iterations. 

Sample comparisons of Navier-Stokes computer 
code results with wind tunnel experimental data are 
shown in figure Bl. In general, these comparisons dem- 
onstrate the capability of the Navier-Stokes method to 
predict the flow characteristics and forces for two- 
dimensional airfoils over a range of free-stream condi- 
tions. As with any computational analysis, the method 
used in this study is subject to some limitations and 
uncertainties. The solution is limited by the numerical 
accuracy of the computer code; CFL3D is accurate in 
time to the second order and in space to the third order. 
Discretization of the geometry results in loss of resolu- 
tion, which can further affect solution accuracy. Proper 
modeling of turbulent flow presents another potential 
problem. The transition point of the boundary layer from 
laminar to turbulent is not known, and how well the 
Spalart-Allmaras turbulence model performs in regions 
of separated flow is uncertain. 

On completion of the limited verification exercise, 
the Navier-Stokes computer code was used to provide 
supplementary information to better define the compli- 
cated dependence of C p Vim on Mach number at low 
speeds. For this purpose the NACA 0009 airfoil was 
evaluated for Mach numbers of 0.05, 0.10, 0.16, 0.20, 
0.30, and 0.50 at a Reynolds number of 9 x 10 6 . The lim- 
iting pressure for the computer code data was defined by 
use of the same iteration process that was applied to the 
wind tunnel data. An example of a plot that was used is 
shown in sketch R. A curve generated by use of these 



Sketch R 


21 



data is superimposed on the plot in figure 7(d), which 
presents data for wind tunnel test results at thfe same 
Reynolds number. 

Navier-Stokes results were also used to provide sup- 
plementary information on the dependence of the limit- 


ing pressure on Reynolds number. For this purpose the 
NACA 0009 airfoil was evaluated for additional Rey- 
nolds numbers of 1 x 10^, 100 X 10*\ and 1000 X 10^ at 
a Mach number of 0.16. Limiting pressure coefficients 
obtained from these data are shown in sketch K. 


O Experimental data (ref. 12) 

Navier-Stokes computer code 



02468 02468 02468 

a, deg a, deg a, deg 


(a) NACA 64-010 airfoil. 

Figure Bl. Experimental and Navier-Stokes computer code (CFL3D) data for two-dimensional airfoils. 
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Figure 1. Attainable thrust prediction method development. 
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Experimental data (ref. 12) 
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a, deg a, deg a, deg 

(a) NACA 64-010 airfoil. 

Figure 6. Example of two-dimensional experimental data to define limiting pressures. 
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O Experimental data (ref. 14) 

— — — Theoretical data, full thrust 
Present method, j- m by iteration 

NACA 0006 


M n = 0.05 M n = 0.11 M n = 0.22 

R n = 1.9 x 10 6 R n = 3.7 x 10 6 R n = 6.8 x 10 6 



0 5 10 0 5 10 0 5 10 

a, deg a, deg a, deg 


(b) NACA 0006 and NACA 0008 airfoils. 
Figure 6. Continued. 
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O Experimental data (ref. 10) 

— — — Theoretical data, full thrust 
Present method, ^ by iteration 

NACA 0006 



8 12 16 0 4 8 12 16 0 4 8 12 16 

deg a, deg a, deg 

(c) NACA 0006, NACA 0009, and NACA 0012 airfoils; M = 0.16. 


Figure 6. Concluded. 



O Experimental data 
Curve fit (eq. (3)) 



(a) R n = 1 x 10 6 . 



(b) R n = 3 x 10 6 . 

Figure 7. Limiting pressure coefficient definition from two-dimensional experimental data. 
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i/c = 0.06 root; 0.04 tip 


O Experimental data (ref. 20) 
Theoretical data without thrust or vortex 
— — — Both methods 

Theoretical data with thrust and vortex 
■ - ■'■■■ "" Previous method 
Present method 






(a) 5 L n = 0°; h Tn = 0°. 

Figure 9. Theoretical and experimental data for 60° swept delta wing fighter with rounded leading-edge NACA 
64A00X airfoil. M = 0.4; R = 2.5x 10 6 . 


41 




i/c = 0.06 root; 0.04 tip 


O Experimental data (ref. 20) 

Theoretical data without thrust or vortex 

Both methods 

Theoretical data with thrust and vortex 
■ - 1 Previous method 

Present method 




Figure 9. Concluded. 
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X/c = 0.06 root; 0.04 tip 


O Experimental data (ref. 20) 

Theoretical data without thrust or vortex 
— — ■ — Both methods 

Theoretical data with thrust and vortex 

■ Previous method 

— — Present method 




(a) h Ln = 0°; h Tn = 0°. 

Figure 12. Theoretical and experimental data for 44° swept trapezoidal wing fighter with NACA 64A00X airfoil. 
M = 0.4; i? = 1.9 x 10 6 . 


45 




O Experimental data (ref. 20) 
Theoretical data without thrust or vortex 

Both methods 

Theoretical data with thrust and vortex 
■ - — — - Previous method 
— - Present method 







Full thrust 


Full thrust 


Figure 13. Theoretical and experimental data for generic arrow wing supersonic transport. M = 0.25; R = 4.8 x 10 ( 







i/c = 0.04 


O Experimental data (ref. 22) 

Theoretical data without thrust or vortex 

Both methods 

Theoretical data with thrust and vortex 
— — - - — — Previous method 

■ " ■ Present method 




(a) Sharp leading edge. 


Figure 14. Theoretical and experimental data for flat modified arrow wing-body combination with various wing sec- 
tions. M = 1.6; R = 2.6 x 10° 
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Experimental data (ref. 23) 


Theoretical data, present method 
With attainable thrust and vortex force 


Without thrust or vortex force 


1 /c = 0.02 root; 0.06 tip 
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O Experimental data (ref. 23) 

Nv 

Theoretical data, present method 

x/c = 0.02 root; 0.06 tip 

Without thrust or vortex force 




Full thrust 


Full thrust 


(c) Large leading-edge radius. 
Figure 15. Concluded. 






O Experimental data (ref. 23) 
Theoretical data, present method 


\ I — — — With attainable thrust and vortex force 

— — — — — Without thrust or vortex force 

x/c = 0.02 root; 0.06 tip 



(a) R = 2.4 x 10 6 . 


Figure 16. Theoretical and experimental data for 61.7° swept delta wing body combination at several additional 
Reynolds numbers. Small leading-edge radius; Af= 0.22. 
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